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LATERAL  STABILITY  OP  AN  AIRCRAF^^AND  VIBRATIONS  OF  THE  AILERONS  WHILE 
TAKING  INTO  CONSIDERATION  ELASTIC  DEFORMATION  OF  THE  WINGS  AND  THE 
ELASTICITY  OP  THE  CONTROL  SYSTEM 


BY  Jerzy  MarYniak,  and  Maria  Ziocka  ((Warsaw) 


1 . Introduction 

In  this  work  there  has  been  investigated  the  influence  of  stiff- 
ness and  damping  in  the  system  of  control  (guidance')  by  ailerons 
while  taking  into  consideration  the  elastic  deformation  of  the  wings 
on  the  lateral  stability  of  the  aircraft  as  well  as  the  vibrations 
of  the  ailerons.  The  aircraft  has  been  treated  as  a rigid  mechanical 
systen^ith  elastically  deformable  wings  and  movable  ailerons. 

The  equations  of  motion  have  been  worked  out  in  quasi-coordinates 
using  the  Boltzmann-Hamel  equations  /i|/  for  mechanica''  'terns  with 
holonomic  constraints  in  a system  of  coordinates  con.  with  the 

aircraft. 


In  the  work  the  assumption  has  been  made  that  the  aerodynamic 
forces  and  moments  do  not  have  any  influence  on  the  form  and  frequen- 
cy of  the  free  vibrations  of  the  wings.  This  assumption  has  made  pos- 
sible a special  examination  of  each  inherent  form  of  wing  vibrations. 
The  v.’ings  constituting  a continuous  system  of  an  infinite  number  of 
dcc'r’oes  of  freedom  have  been  substituted  for  by  a precisely  determined 
num  )or  of  dotrees  corres'-'onding  to  the  number  of  assumed  forms.  The 
'’nr  and  frequencies  the  free  ■'Mbrations  have  been  determined 
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experimentally  by  way  of  resonance  investigations ./6, 8/. 

The  lineari 7ation  of  the  equations  of  motion  v;as  carried  out  on 

the  basis  of  the  theory  of  small  perturbations  /I ,2 , 3 , 7 ,20/.  It  was 

( 

assumed  that  antisymmetric  motions  of  the  aircraft  cause  exclusively 

antisymmetric  changes  of  the  aerodynamic  forces  and  moments,  moreover 

symmetric symmetric  changes  of  the  aerodynamic  loads.  The  above 

assumptions  made  possible  the  resolution  of  the  system  of  equations 

/I ,2,3,204  which  describe  an  arbitrary  motion  of  the  aircraft,  into 

two  systems:  the  system  of  equations  of  symmetric  motions  /8,9,13/ 

antisymmetric 

of  longitudinal  stability,  and  the  system  of  equations  of  motions 

\ 

of  lateral  stability 

There  have  been  taken  into  account  five  degrees  of  freedom*,  of 
these  three  are  the  degrees  of  freedom  of  the  rigid  aircraft:  roll'll 
yaw  the  velocity  of  lateral  di  snlacements  v as  well  as  the  anti- 
symmetric elastic  deformations  of  the  wings  C and  the  elastic  deflec- 
tion of  the  ailerons 

After  linearization  of  the  system  of  equati ons, the  solution  was 
carried  out  as  far  as  the  designating  of  the  eigenvectors  and  the 
eigenvalues  of  the  matrix  of  state  corresponding  to  them.  Sample  nu- 
merical computations  have  been  carried  out  for  an  airplane  of  the 
';vilgs"  class  in  accordance  with  their  own  programs  in  the  Establish- 
ment for  Numerical  Calculations  of  the  University  of  Warsaw. 

In  the  accessible  literature  having  to  do  with  the  dynamics  of 
moving  object s , one  does  not  meet  with  a derivation  of  the  eiuations 
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of  inotlon  thru  use  of  the  Hamel-Bolzmann  equations  for  systems  of 
holonomic  constraints.  Properly^ the  use  of  the  Boltzmann-Hamel  equa- 
tions for  the  deriving  of  the  equations  of  motion  of  mobile 
(moving)  objects  in  a system  of  coordinates  connected  to  that  object 
makes  possible  in  a relatively  simple  manner  a taking  into  conside- 
ration of  the  degrees  of  freedom  that  result  from  the  deformabi li ty 
of  a body  as  also  the  relative  motions  of  the  elements  of  the  object 
being  examined 

The  Bo Izmann -Hamel  equations  are  generalized  Lagrange  equations 
of  the  second  type  expessed  in  quasi-coordinates  and  quas i -veloci ties. 
The  quasi-velocities  are  linear  combinations  of  the  generalized  velo- 
cities whose  cooeff icients  are  dependent  on  the  generalized  coordi- 
nates /l4.,19/,  where  one  can  include  also  free  (arbitrary)  terms  as 
also  clearly  (those)  dependent  on  time.  /19/.  In  the  case  considered 
such  quasi-velocities  are  kinematic  parameters  of  motion  determined 
in  a sysfem  of  central  axes  of  reference, rigidly  connected  with  the 
aircraft.  The  kinematic  parameters  mentioned  above  are  the  angular 
velocities  of  the  aircraft  P,Q,R  as  well  as  the  linear  velocities 
of  its  center  of  mass  U,V,W  /I  , ? , 3 ,4*'^ » 3 , 1 )4 , 1 9 , and  20/. 

The  equations  of  motion  introduced  in  the  third  section  of  the 
present  work  are  universal  and  one  can  anply  them  directly  to  the 
deacrlntion  of  the  motion  of  arbitrary  deformab?§^iA5®cts  in  the  as — 
sumod  system  of  reference. 
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2 .The  Assumed  System  of  Reference 


For  the  description  of  the  dynamics  of  an  aircraft  three  systems 
of  reference  are  Indispensable:  the  gravitational  system  precisely 
connected  witi<  the  earth  Ox-|y^z-|,  the  velocity  system  connected  with 


the  flow  Ox  y ? as  well  as  the  one  connected  rigidly  with  the  air 
a a a^ 


craft  Oxyz . 


Diag  #1  The  assumed  system  of  reference  Oxyz  connected  with  the  air 
craft  as  well  as  the  introduced  linear  and  angular  velocities. 


The  momentary  situation  of  the  aircraft  as  a rigid  body  is  de- 
termined by  the  situation  by  the  center  of  mass  of  the  ob.iect  r,(r. . , -.i 

measured  with  respect  to  the  immovable  system  of  coordinates  Ox^y^z^ 
connected  with  the  earth  as  well  as  of  the  angles  of  rotation  of  the 
aircraft  , 0, 


The  angles  of  rotation  determine  uniquely  the  situation  of  the 
system  of  coordinates  precisely  connected  with  the  aircraft  Oxyz  with 


respect  to  the  gravitational  svstem  of  coordinates  Ox  v z^which  is 

F'  g R 


parallel  to  the  immovable  system  Ox^  y^  ( (diag  ) . 


. 
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The  assumed  angles  of  rotation  are  quasi-Euler  angles,  called 
also  aircraft  (angles )/1 ,2,1  9/. 

The  names  of  these  angles  are  as  follows;  </>  the  angle  of  roll, 

0 -the  angle  of  pitch,  '/■'  -the  angle  of  yaw. 

The  motion  of  an  aircraft  was  described  in  the  central  system 
0xy7.  rigidly  connected  with  the  aircraft,  around  axes  directed  as  on 
diag  #1  and  diag  #2. 


*9  'H'L 

'T' 


Diai^ram  #2  ^Assumed  components  of  forces  and  momenta  in  the  reference 
system  Oxyz  connected  with  the  aircraft. 


The  components  of  the  momentary  velocity  vectors,  linear  i-'  and 

ancrul  ar /■?  in  the  assumed  coordinate  system  ((diag  #1)  are  as  follows: 

vector  of  linear  velocity 

(1  ) i'c  tVl  Vj\  U’k, 

where  U doaignates  the  longitudinal  velocity,  V the  lateral  velocity 
and  W--the  velocity  of  vertical  (plumb)  displacements. 

--the  vector  of  momentary  angular  velocity  rj 


5 


where  P is  the  angular  velocity  of  roll,  Q--tho  angular  velocity  of 
pitch,  R--the  angular  velocity  of  yaw. 


The  vectors  of  the  external  forces  and  the  moments  of  external 
forces  iscttng  ^he  aircraft  have  the  form  (diag  #2): 

vector  of  the  external  forces  y 

(3)  F 

where  X signifies  the  longitudinal  force,  Y-«-the  lateral  force,  7.-- 
the  vertical  force 

the  vector  of  the  main  (overall)  moment  j( 

(I4)  .rr  ^ i/\  .\fj  > Xk", 

where  L is  the  roll  moment,  M--the  pitch  moment,  and  N--the  yaw  mo- 
ment . 

The  angular  velocities  P,Q,R  are  linear  combinations  of  the  gene 
ralized  velocities  and  V'  with  coefficients  dependent  upon  the 
generalised  coordinates  </>,  0 and  V'  and  are  expressed  in  the  followitig 
form : 


F 

1 

0 sint'l 

(^) 

f) 

0 

1, s!,iO>  f* 

n 

0 

'j 

' c 

0 

■ '\.  ^n 

V' 

»v 

1 

The  kinenat'c  relations  between  the  linear  velocities 
measured  in  the  non-moving  system  Ox^  y^  7 and  the  components  of  velo  — 
cityU,V,W  are  the  following: 


The  relations  (1)  and  (6)  determine  the  kinematic  nararneters 
that  are  quasi-velocities. 


3 . Equa  " Motion  of  a Deformable  Free  Object 

The  model  of  a non-deformable  aircraft,  most  frequently  encoun- 
tered in  the  litterature,  cannot  always  be  assumed  (accepted)  in  an 
;nvestip;0t^ on  of  the  dynamic  properties  of  an  ob,)ect.  Some  T«ip;id  and 
elastic  T«elative  motions  can  have  a substantial  influence  on  the 
'.haracter  of  the  motion  of  an  aircrart.  In  the  case  of  taking  into 
consideration  the  elastic  flexibility  of  the  wings, one  obtains  a sys- 
tem of  an  infinite  number  of  degrees  of  freedom  of  motion.  Practical 
carrying  out  (if  the  computations  for  such  a system  is  impossible  and, 
therefore^ there  is  also  used  an  anoroxlnate  method.  It  1s  based  on 
the  assumption  that  the  aerodynamic  forces  and  moments  do  not  change 

1 

the  non-elastic  free  vibrations  of  the  wings. 

[ 

The  sag  of  a wing  (dlag  7/ii ) is  described  by  the  function:  ( 

^...0  j 

II 
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r 


Dlag  #3:  Assunied  model  of  the  aircraftand  of  the  angular  displace- 

ments of  the  antisymmetric  rotations  of  the  aircraft. 


r>-!.  4.  nuiJcl  nichow  pr,-ccliyl,ij:;cych  i oJ'K.’taK'ai  rivtnydi  sVr,-y4cl  sainolotu 

Diag  Assumed  model  of  the  rolling  motions  and  elastic  deforma- 

tions of  the  wings  of  the  aircraft. 

where  h.(y)  is  the  successive  form  of  the  free  vibrations.  This  has 
made  possible  an  examination  of  the  inf luence  on  the  motion  of  the 
aircraft  «f each  form  of  the  vibrations  separately.  In  agreement  (con- 
formity) with  the  above  that  correspond  to  the  ith  form  have  been  pre- 
sented as  follows: 

(7)  'Vy.O  ’:S)V-X0- 
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There  has  been  taken  into  consideration  also  the  motion  of  the 
ailerons,  which  is  possible  in  spite  of  the  blocked  control  stick, 
due  to  the  existence  of  elastic  deformations  in  the  system  of  control 
by  ailerons.  The  displacement  of  the  ailerons  is  determined  by  the 
angle  of  rotation  of  the  aileron  ground  the  axis  of  the  hinges(diag  3) 


The  equation  of  motion  of  the  aircraft  has  been  introduced  in 
quasi-coordinates,  using  the  Boltzmann-Hamel  equations  for  holonomic 
systems  /J4./.  The  Boltzmann -Hamel  equations  are  generalized  Lagrange 
equations  of  the  second  type  for  non-inertial  systems  described  in 
quasi-coordinates,  and  they  have  the  following  form: 


where  1.1.  . ,/,,A  signifies  the  number  of  degrees  of 

freedom,  --the  quasi-velocities,  T-"---the  kinetic  energy  in  quasi- 
velocities, "^^--the  quasi -coordinates  , --the  generalized  forces. 


The  relations  between  the  quas i-veloci ties  and  the  generalized 
velocities  has  the  form: 


(Q) 

(10) 


k 

ff  I 

k 


- 


/'♦ 


where  ,'r^  signifies  the  generalized  velocities,  ,r.,  g..  , 

the  generalized  coordinates,  , g:,  ...gO,  whereby 

there  exists  the  follov;ing  matrix  dependence  * 


(11  ) 


l^,.l  '■ 
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The  triple  indexed  multipliers  (coefficients)  of  Boltzmann 
are  determined  by  the  relationship: 


k k 


(12: 


■ >’  >'('>  ■ ■ 
<■  ‘^1  < I ’ 


1 I j . 1 

In  the  case  when  the  quasi-coordinates  are  generalized  coordinates, 
then  the  triple- indexed  Boltzmann  coefficients  7'..  (12)  ai  e equal 
to  zero. 


In  the  assumed  model  of  a deformable  aircraft  the  vector  of  the 
quasi-velocities  is  the  following: 

where  U,V,W,P,Q,R  determine  the  relations  ( ')  and  (6),  and  the  vector 
of  quasi -coordinates  corresponding  to  it  has  the  form: 

( 1 ) <£  -=  ool[.rt-,  .T,-.  :rp,  :tq,  .r.,;,  /?,  c]. 

The  vecibr  of  the  generalized  coordinates  is  as  follows: 


(15) 


<1  /?.n- 


The  matrix  [•,  .]  in  the  case  of  the  assumed  model  (diag#3  and 
diag  #14.)  in  the  assumed  system  of  coordinates  (diagvl  ) in  conformance 
with  the  relations  (5)  and  (6)  has  ^he  following  form: 


Ay  \ 

0 ! 

0 

(16) 

[r,,,]  = 

0 

: 

0 

0 ; 

0 ; 

I 

whereby  the  matrix 

corresponding 

to  it 

(17) 


Ap'  0 0 

0 0 

0 0 I 
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The  majority  of  the  Boltzmann  coefficients  y',_  in  the  case  of 
the  assumed  deformable  aircraft  are  equal  to  zero.  One  of  the  non- 


zero Boltzmann  coefficients  is  yU-  It  is  computed  in  accordance  with 
the  following  relationship: 


8 » , \ 

\ ’ \ 1 / (>6I  \ , . 

= > > I . - , 

7,  <.7.  I 


(18) 


'A* 


I 1 


The  quasi -coordinates  & C are  generalized  coordinates.  In  conformity 
vjith  the  above 


0;  0 


for  dj  1,2, 8,  is  simplified  and  the  relationship  (I8)  assumes 

the  form 


a -■  I ^ 1 


/^\cqi  <-<h  I 


Analyzing  the  matrices  ' and  rermits  us  to  note  that: 

6|4  = /?24  ~ ^34  “ ^54  ~ ^64  ~ ^44  1* 

<^61  “ ^f>2  ~ *^63  ” ^64  “ 0;  fti5  = 6j5  = bjf  = 0, 


stqd 


/ie 


1^4 


‘/s*  --  - f'fsf’i*-  ~ c<,s'’-:4. 

C(j4 


d(  Jiny) 

- ■ = --  cOby 

dq4  d(p 

Scet  (?(cosy  cos(9) 


pi  /.y  czyin 


dq4 


C(p 


= — siny  cObO, 


W rc/iiltacic  otr/.jniano 


his  - /»6 


sinyi 


cob(y 


yf4  cos’v  I sin^y  - I . 


The  regaining  non-zero  Boltzmann  coefficients  are  computed  anaiogous- 

y\b  — ~ I.  >'(.3  li 

V*  1 


ly:  they  have  the  following  values: 


t 


Af'.cr  introducing  the  thus  computed  Boltzmann  coefficients  into  the 
equation  (8)we  get  equations  of  motion  for  an  arbitrary  free  obiect 
whose  motion  is  described  in  the  assumed  system  of  reference.  The 
equations  of  motion  take  on  the  following  form: 


(.19) 

d lor*\ 

or* 

or* 

or* 

/ 

dt  \ cu  ) 

O.iu 

cy  ^ '• 

c^yQ 

-Qt. 

(^0) 

dt  \ C K / 

or* 

Orty 

1 

‘I'  * />  - 

t Ik 

d lor*\ 

or* 

cT*  ^ 

or* 

(21) 

~dt  \ ay) 

cU 

oy^^ 

= Qi; 

(22) 

'n*\ 

1 or* 

or* 

or*  „ 

or* 

or* 

,1,  I 

O.Tp 

— 

oV 

y 

av 

• 

OR  ^ 

(23) 

d 1 

'cn\ 

6T* 

or* 

or* 

or* 

( 

1 dr* 

'll/  ^+ 

olV 

OP 

p - 
OR 

’(?5, 

(21) 

d , 

IdT*' 

( dT* 

or*  „ 

or*  „ 

or*  ' 

or* 

~dt  * 

I 

- cu 

p 

cQ 

(25) 

d 

'dt 

1%)- 

or* 

0?' 

Oh 

(26) 

d 

ior*\ 

or* 

dt 

( ot  ] 

-Qt- 

The  equations  (19)  & (2ii)  describe  the  motion  of  an  arbitrary 
rigid  body  in  a central  system  of  coordinates  connected  with  the  ob- 
ject. The  remaining  two  equations  are  the  result  of  taking  into  con- 
sideration additional  degrees  of  freedom;  of  the  relative  motions  of 
the  ailerons  (2'5)  and  the  deformation  of  the  wings  (26). 

In  an  arbitrary  motion  of  an  arbitrary  deformable  object,  the 
number  of  equations  of  the  type  (26)  can  he  arhitrslV  and  depends  ex - 
clusively  on  the  number  of  additionally  taken  into  consideration  de- 
grees of  freedom  in  the  presence  of  the  unchanging  form  of  the  first 
six  equations  (19)-(2l|). 
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It.  Squatirna  of  the  Antisjmnetro  c!  Motions  of  a Deformable  Aircraft 

In  an  arbitrary  motion  of  an  obieot,  equations  (19)  thru  (26)  in  ge- 
neral do  not  resolve  themselves  into  equations  describing  symmetric 
and  antisymmetric  motions.  They  are  very  non-linear  ordinary  differen- 
tial equations  of  the  second  order.  The  resolving  of  the  equations 
is  only  possible  in  use  for  studies  of  the  theory  of  small  perturba- 
tions with  regard  to  an  established  motion  end  linearization  of  the 
equations . 


In  the  present  work  we  have  assumed  that  the  aircraft  carries 
out  only  antisymmetric  motions  (diag  3 and  diag  #Ij. ),  that  means  the 
yawing  motion  the  rolling  motion  the  lateral  displacement 
the  antisyrimetric  displacements  of  the  ailerons  p and  the  antisymme- 
Trio  vibrations  of  the  wings  ^ In  the  presence  of  the  above  assunp- 


(27) 

(28) 
(29) 
(’0) 
(31) 


system 

of  five  equations 

in  the 

general 

d ICT*\  CT* 

CT* 

CT* 

'dt  \ C'  F / e.rty  + 

cu^~ 

d lvT*\ 

or*  CT*  CT*  „ 

CT* 

CT* 

di  \ C P 1 

(-T,  cV  c\V 

,-a 

d /cr*\ 

1 CT*  CT*  „ CT*  „ 

CT*  ^ 

^T* 

dl  \ cR  1 

CTlfi  CU  OY 

CP 

igP‘01. 

CT* 

dt  1 

1 ‘V  ) 

1 cp 

d 

, CT* 

dt 

1 rC  ] 

1 cc 

0}, 

Qt- 


The  whole 
velocities  has 


kinetic  energv  of  the  aircraft  T-t-  computed  in  quasi.— 
the  followin'*  form: 


1 


(32) 


r*  = 2*  H KM  >3'*)  + /,  C'  f /.  I'^  * /, + /.  K']  - M ^’2  <- 

+ S,{iVP-UR)  i S,iVR~  UQ)  i SXUQ-yP)  I ((-><3  t B'  \ B'i)\V+ 
+ {A,  + b:  yB'i)P-[A^  HB1  i 7?I).v^  ~b:-B',]Q}C  \ 

where 

j bJ2  ‘/j 

-kl2  -0/2 

fr/J  <’/2 

Ai  f ^4  = / •‘^/.'■)''(3)'(>’. 

- */2  l>,'2 

*,'2  fc/Z 

i/2  it  *.-  *■- 

i/2  i/2 

i/2  it  '■.'2  it 


where  MgSigntries  the  mass  oT  the  whole  aircraft,  l J.J  ore 

moments  of  inertia  and  the  moment  of  the  deviation  of  the  aircraft, 
with  respect  to  the  system  of  reference,  Oxyz : S ,S  ,S  --are  the  static 

moments  of  the  aircraft  x^ith  respect  to  the  reference  system  Oxyz: 

R,  f’,  B,,  y,,  Si  --are  the  moments  of  inertia, deviation  and  static  of 

the  ailerons  with  respect  to  the  axis  of  the  hinges  v and  of  the  axis 
of  symmetry  of  the  aircraft  f,  ; the  upper  indices  of  h & P determinf 
respectively  the  left  and  right  aileron;  are  the  distri- 

butions of  the  masses  of  the  wings  and  of  the  aileron  as  a function 
Span;  h(y)  — is  the  function  of  the  sag  of  the  wing  corresrondi nr 
to  the  form  being  investigated  of  the  free  vibrations. 

The  generalized  forces  anpearing  on  'he  right  sides  of  the  euna- 
tions  i?-7)  thru  ( 31  ) determine,  taking  into  consideration  the  r oten— 
tial  enerry  of  the  deformations  of  the  winrr-  and  of  the  eon'.rol  sys' 


F 


system  [8,  the  gravitational  forces  £l,  2, 

as  the  aerodynamic  forces  and  moments  [1,  2,  3, 


3,  8,  as  well 

8,  14,  20J. 


The  potential  energy  of  the  elastic  deformations  of  the  wings 
and  of  the  control  system  of  the  ailerons  has  the  following  form: 


(33) 


U. 


Where  the  elastic  rigidity  of  the  wings,  the  stiffness  of 
the  system  of  control  of  ailerons,  where: 


(34) 


6/2 


/ '’h{y)k-{y)^>y 

k/2 


is  the  generalized  stiffness  of  the  wings  that  corresponds  to 
the  taken-lnto-conslderatlon  form  of  the  free  vlbratioiis,  which 
(stiffness)  is  described  by  the  function  h(y)  of  frequency  of 
vibrations  oj. 

The  viscous  damping  in  the  system  of  control  by  ailerons 
has  been  taken  into  account  by  the  Introduction  of  Rayleigh's 
dissipation  function  U„. 

n 


(35) 


The  components  of  the  force  of  gravity  in  the  reference 
system  Oxyz  have  the  form 


(36) 


• '>‘S 


where 


A. 


-sin  6) 
cosOsin'.f*  , 
cos  O I,  Ob  0 
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and  In  the  case  under  consideration,  which  takes  into  consideration 
exclusively  anti-aymmetric  motions 

^37)  y,  — «ji’cosf9siii0. 

The  aerodynamic  forces  and  moments  acting  on  the  aircraft  have 
been  Introduced  while  taking  into  consideration  the  stationary 
aerodynamics.  Linearization  of  the  aerodynamic  forces  and  moments 
has  been  carried  out  in  accordance  with  the  method  of  Bryan  [1,  2, 

3,  20].  This  method  is  based  on  assumption  that  the  aerodynamic 
moments  and  forces  are  momentary  functions  of  the  magnitude  of  the 
changes  of  velocity,  linear  and  angular  and  of  their  derivatives. 
These  functions  have  been  developed  into  a Taylor  series  with  respect 
to  the  previously  mentioned  changed.  In  these  series  there  have 
been  considered  only  terms  of  the  first  order  [1,  2,  3,  20]. 

In  the  case  examined  it  has  been  assumed  that  the  aircraft 
moves  with  a steady  uniform  level  motion.  One  assumes  that  the 
steady  motion  of  the  aircraft  is  subjected  to  small  perturbations, 
that  means  that: 


0 . 

<r> 

P 

P> 

U 

■-  Uo  const, 

0 . 

- const, 

Q 

0, 

K 

>J' 

^ V. 

R 

r. 

U' 

0. 

The  generalized  forces  in  the  equations  (27)-(31)  when  taking 
into  consideration  the  above  perturbations  and  the  introduction  of 
the  relations  (37)  and  (38)  have  the  form: 

Qy  ■ ^ ! fpP  ) K, r 1 

Qf  t aV  I fp/J  1 J I /;C  I 

( 39 ) Q%  'V.-  I i .V,^  ! .V;/?  I ,V;C  1 A',/?. 

Of  RpP  \ \ R;,;i  ! /?,•;  i 

Q*  F,P  ! p,r  \ F;ji  ! /•;.•;  ! 

where  y Y --  ‘ E-  ■ ‘ according  to  the  names  accepted  In 

' 2v  ' " Pp  ' ' c' 

flight  lore,  the  aerodynamic  derivatives 

Id 


[1 


2,  3,  20]. 


The  aerodynamic  derivatives  that  appear  in  the  relations  (38) 
were  Introduced  in  [1,  2,  3,  6,  8,  1^,  and  20]. 

The  system  of  equations  (27)  thru  (31),  after  taking  into 
consideration  (28)  and  (39),  are  converted  into  a non-dimensional 
form  dividing  the  equations  of  the  forces  by  oi'-S,  and  the  equations 
of  the  moments  by 


cv^s 


•x 

AUbizy 


I - 


t 

A 

t 


V 


P = pi 


y,  =-■ 


c.sr,ft/2 


N, 

C.S-K,A/2 


fr 


QS\\hl2 


the  aerodynamic  time, 

the  relative  density  of  the  aircraft, 
the  non-dimensional  moment  of  inertia, 

the  non-dimensional  time. 


non-dimensional  linear  velocity, 

non-dimensional  angular  velocity, 

non-dimensional  derivative  of  lateral  force  with 
respect  to  the  change  of  the  linear  velocity  of 
sideslip , 

non-dimensional  derivative  of  the  lateral  force 
with  respect  to  the  change  of  roll  angular 
velocity , 

non-dimensional  derivative  of  the  yawing  moment 
with  respect  to  the  change  of  the  linear  velocity 
of  sideslip, 

non-dimensional  derivative  of  the  roll  moment  with 
respect  to  the  yaw  angular  velocity  change. 


In  an  analogous  manner  there  have  been  presented  in  a non- 
dimensional  form  the  remaining  terms  (expressions)  of  the  equations. 


The  system  of  equations  of  motion  in  non-dimensional  form  has 
been  obtained  in  matrix  notation  in  the  following  form: 
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where 


Ax  1 I>x  1 Cx  0, 

X - = col[.Ti.,  rr,,  P], 

X = r,  C,  /*], 


whereby 

— A is  the  matrix  of  the  coefficients  of  Inertia,  VIZ 


1 0 0 0 0 

0 1 l-Jj. 

0 -jJh  1 0 0. 

0 0 1 

-0  Ji-Jh  0 /•{•/y,  1 


— B is  the  matrix  of  coefficients  of  damping,  VIZ 


-y. 

-ypift 

1 0 

0 , 

-I'ljjx 

-IplJx 

-4/A 

^-/;/A 

-IplJ. 

-IrlJz 

-"•/A 

-npUz 

0 

^p 

fr 

.‘’/i 

0 

~ ’^pOi 

-'rU,, 

-^ik 

Jn 

— C is  the  matrix  of  coefficients  of  stiffness,  VIZ 


0 -y^  0 0 0 

0 0 0 0 -V//, 

0 0 0 0 

0 0 0 Cfi 

0 0 0 0 


The  matrix  equation  (^0)  of  second  order  gets  reduced  to  an 
equation  of  the  first  order  in  the  form: 

(41)  i’«i : Q'l  0.  q I ], 
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where  the  matrix  of  state  R has  the  form 


« I .f  I'll 


The  solution  has  been  reduced  to  the  designation  (evaluation) 
of  the  eigenvectors  and  of  the  eigenvalues  corresponding  to 

‘-‘'ij.jn  of  the  matrix  of  state  R (^3) 

The  general  solution  has  the  form 

I 

( 44  ) *1(0  = 2j 


where  the  Cj  designate  constants  dependent  upon  the  initial  condi- 
tions which  are  the  values  of  the  perturbations  from  steady  motion 
for  the  moment  t=0,  hj  is  the  frequency  of  oscillation  of  period 

T t,  . is  the  coefficient  of  damping,  whereby  ^ P 


is  the  time  of  damping  down  of  the  amplitude  to  ^ for  f < 0,  and  in 
case  I > 0,  the  time  of  doubling  of  the  amplitude. 

5.  Numerical  Example  and  Suggestions 


The  example  computations  have  been  carried  out  for  a light  air 
craft  of  the  tourist  class  "Wilga".  We  resolve  the  system  of 
equations  (40)  evaluating  the  eigenvectors  q^j  and  the  corresponding 
eigenvalues  Xj  of  the  matrix  of  state  R (43). 

All  the  computations  were  carried  out  in  accordance  with  their 
own  programs  at  EMC  GIER  in  the  Institute  of  Numerical  Calculations 
of  the  University  of  Warsaw. 
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— Figui-t;  5.  Coefficients  of  damping 
^ and  frequencies  of  oscillation 
n of  the  first  four  eigenvalues 
X-,  through  Xji  as  a function  of 
■^the  stiffness  of  the  system 
of  control  of  ailerons. 


The  results  have  been 
presented  in  the  form  of  graphs 
(Figures  5-9)  on  which  there  have 
been  represented  by  the  continuous 
line  changes  of  the  damping  co- 
efficients Cj  and  with  a broken 
line  the  frequency  of  oscillation 
rij  . There  have  been  designated 
with  the  same  Indices  on  all 


graphs  the  eigenvalues  correspond- 


ing to  them, /that  characterize  the  same  motions  of  the  aircraft. 


wings  and  ailerons: 


t 

s 1 


^■2 


^2 


X.3,4  ■ ■ s J,4  1 


or 

h - f* 

^■7S  ' S7»  I <’/7,S 


characterizes  aperiodic  displacements  of 
the  ailerons  6 always  powerfully  damped 

Ci<0, 

characterizes  the  aperiodic  spiral  motions 
that  indicate  weak  instability 

characterizes  periodic  oscillations  that 
correspond  to  rocking  turns  p and  v coupled 
with  a yawing  motion  r of  a motion  always 
damped  4^*^ 

characterizes  motions  periodic  or  aperiodic 
of  the  ailerons  6 coupled  with  the  rolling 
motions  of  the  aircraft  p 

damped  ^^g<0,or  divergent  ^>0 

characterizes  elastic  vibrations  of  the  wings 
always  damped  ^yg<0  of  frequency  n^g  close 

to  the  frequency  of  the  free  vibrations  of 
the  wings  u. 


I 
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a)  Influence  of  the  stiffness  of  the  control  system  on  the 
stability  of  the  airplane.  An  Increase  of  stiffness  in  the  control 
system  (Figures  5 and  6)  causes  a decrease  in  damping  of  aperiodic 
displacements  of  the  ailerons  (3)C2  with  unchanged  characteristics 
of  the  spiral  motions  and  of  the  rocking  turn  (Figure  5). 

An  increase  in  the  stiffness  also  does  not  influence  the  elastic 
vibrations  of  the  wings  nor  the  frequency  ri.^g  as  also  the  damping 
^yg<0,  it  has,  however,  a powerful  influence  on  the  damping  and 
character  of  the  motion  of  the  ailerons  that  is  coupled  with  the 
roll  of  the  aircraft  (Figure  6). 


Figure  6,  Coefficients  of 
damping  ^ and  of  frequency 
of  oscillation  q of  the 
eigenvalues  Ag-Xg  as  a function 

of  the  rigidity  of  the  system 
of  control  of  ailerons. 
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With  small  stiffness  of  the  control  system  of  ailerons,  the 
displacements  of  the  ailerons  g and  the  rolling  motions  of  the 
aircraft  p are  aperiodic  diverging  motions  ?g>0  and  5g>0,  which  in 
the  presence  of  an  increase  of  stiffness  pass  over  into  periodic 
motions  of  frequency  initially  divergent  C^g<0,  and  then 

damped  (Figure  6). 

b)  Influence  of  Damping  in  the  System  of  Control  of  Ailerons 
on  the  Stability  of  the  Aircraft.  An  increase  in  the  viscous  damp- 
ing in  the  system  of  control  of  ailerons  (Figures  7 and  8)  causes 
an  increase  of  the  powerful  damping  ^^^<<0  of  the  aperiodic  dis- 
placements of  the  ailerons  g (Figure  7)  with  an  unchanged  character 
of  the  spiral  motions  of  the  aircraft  well  as  of  the 

frequency  the  damping  of  the  rocking  turns  of  the 

aircraft  (Figure  7).  Damping  in  the  control  system  also  does  not 
affect  the  frequency  n,^g  and  the  damping  Cy0<O  of  the  elastic 
vibrations  of  the  wings  (Figure  8). 

A change  in  the  damping  in  the  control  system  has  a decisive 
and  most  important  effect  on  the  displacement  of  the  ailerons  6 
and  coupled  with  them  on  the  rolling  motions  of  the  aircraft  p 
(Figure  8).  In  the  presence  of  a small  damping  there  appear 
harmonic  oscillations  of  frequency  initially  powerfully 

divergent  passing  over  into  damped  with  a simultaneous 

decrease  in  the  frequency  of  oscillation  At  a certain  critical 

damping,  periodic  vibration  of  the  ailerons  and  wings  pass  over 
into  powerfully  damped  aperiodic  motions  ^^<0  and  Cg<0  (Figure  8). 
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B’lgure  ?.  Figure  8. 

Figure  7.  Damping  coefficients  C and  frequencies  of  oscillation 
n of  the  first  four  eigenvalues  as  a function  of  the  damping 

of  the  system  of  control  of  ailerons. 

Figure  8.  Coefficients  of  damping  5 and  frequencies  of  oscillation 
n of  the  eigenvalues  as  a function  of  the  damping  of  the 

system  of  control  with  ailerons. 

c)  Influence  of  the  balance  (trim)  of  the  Ailerons  on  the 
Stability  of  the  aircraft:  Previous  balancing  of  the  ailerons 

effectively  Influences  the  dynamic  properties  of  the  aircraft 
as  also  even  the  motions  of  the  ailerons  themselves  by  the 
stabilizing  of  the  aircraft  (Figure  9). 
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A chanf^e  in  the  balance  has  a decisive  influence  on  the  displace- 
njent  of  the  ailerons  and  on  the  rolling  motions  of  the  aircraft  p. 
Zero  static  balancing  .lust  as  before  (that  means  the  center  of  mass 
of  the  aileroriis  located  on  the  axis  of  rotation  of  the  aileron  or  in 
front  of  the  axis  causes  aperiodic  motions  powerfully  damped 
and'^^^C  , moreover,  rear  balancing  favors  the  appearance  of  periodic 
oscillations  of  frequency /V|^^  and  of  damning 

6,  General  Suggestions 

Taking  into  consideration  of  additional  degrees  of  freedom  which 
are  the  elastic  deformability  of  the  wings  and  the  elastic  displace- 
ments of  the  ailerons^  in  relation  to  the  results  obtained  in  the  case 
of  a rigid  plane,  causes  the  appearance  of  four  additional  eigenvalues. 

The  eigenvalues  precise  equivalents  of  the  eigen- 

values that  characterize  the  motions  of  a rigid  aircraft,  that  is^ 
spiral  motions  and  rocking  turns. 

In  the  case  considered,  there  is  a lack  of  an  equivalent  that  characterizes 
the  aperiodic  powerfully  damped  roll  of  the  rigid  aircraft. 

There  arises,  moreover,  a powerful  coupling  of  the  relative  displace" 
ments  of  the  ailerons*^  with  the  rolling  motions  of  the  aircraft  p. 

Pov;erfully  damped  elastic  vibrations  of  the  wings  do  not  influence 


in  a substantial  manner  the  remaining  motions  of  the  aircraft; 
they,  moreover,  exclusively  depend  on  the  rigidity  of  the  wings 


The  results  obtained  and  the  suggestions  extracted  on  the 
basis  of  them  are  fair  (proper)  for  the  case  considered.  The 
application  of  them  to  another  type  of  aircraft  or  flying 
object  requires  additional  numerical  computations  in  accordance 
with  the  programs  worked  out. 


Figure  9*  Coefficients  of  the 
damping  C and  of  the  frequency 
of  oscillation  n as  a function 
of  the  level  of  static  balance 
of  the  ailerons. 
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